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We study the properties of hypernuclei containing one A hyperon in the framework of 
the correlated basis function theory with Jastrow correlations. Fermi hypernetted chain 
integral equations are derived and used to evaluate energies and one-body densities of 
A-hypernuclei having a doubly closed shell nucleonic core in the jj coupling scheme, from 
Carbon to Lead. We also study hypernuclei having the least bound neutron substituted 
, by the A particle. The semi-realistic Afnan and Tang nucleon-nucleon potential and Bod- 

mer and Usmani A-nucleon potential are adopted. The effect of many-body forces are 
considered by means either of a three body A-nucleon-nucleon potential of the Argonne 
type or of a density dependent modification of the A-nucleon interaction, fitted to repro- 
duce the A binding energy in nuclear matter. While Jastrow correlations underestimate 
the attractive contribution of the three body A interaction, the density dependent po- 
tential provides a good description of the A binding energies over all the nuclear masses 
range, in spite of the relative simplicity of the model. 
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I. INTRODUCTION 



The microscopic study of atomic nuclei has gained in the recent years a noticeable level of reliability, 
following the improved knowledge of the interaction between the nucleons and the development of new 
and more sophisticated many-body theories. It is now possible to exactly solve the Schrodinger equation 
for nuclei containing up to A = 8 nucleons 0] , while the structure of medium and heavy nuclei can be to 
a large extent understood by means either of variational techniques |^] or of highly refined perturbative 
expansions These relatively recent developments have clearly stressed several inadequacies of the 
independent particle model (IPM), on which the time honoured shell model is based. For instance, the 
; ^ ' IPM completely fails in describing the large energy behaviour of the one nucleon momentum distribution 
I and the quenching of the spectroscopic factors Q . 

The scenario is not equally bright for hypernuclei, or bound systems of nucleons and one or more 
strange baryons, as the A or E hyperons. The many-body techniques developed for the nuclei can be, 
in most of the cases, straightforwardly extended to the hypernuclei, at least from a formal point of view. 
However, the weak point of a consistent microscopic approach now lies in the poorer understanding of 
the hyperon-nucleon (YN) and hyperon-hyperon (YY) interactions with respect to the nuclear case. The 
limited amount of information on the YN scattering data is a consequence of the experimental difficulties 
due to the short lifetime of the hyperons and to the low intensity of the beam fluxes. For this reason 
phenomenological models and effective YN interactions have been often used in the study of hypernuclei. 
Approaches based upon Skyrme interactions Q and Woods-Saxon ||] potentials as well as relativistic 
mean field theories |^-^ have been used. However, also microscopic YN potentials have been recently 
produced [^-^ and adopted to study the structure of hypernuclei from to 2°^Pb ||l3|,U. The A 
single-particle properties have been investigated in details in Ref. ]l5[ | , making use of the nuclear matter 
Fermi hypernetted chain results in a local density approximation. 

In this paper we present a microscopic study of the structure of various hypernuclei containing a 
single A hyperon. This study has been done by extending the correlated basis functions theory and the 
Fermi hypernetted chain (FHNC) equations technique m& to describe nucleonic systems containing one 
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hyperonic impurity. The nucleon-nucleon (NN) interaction we have chosen is the semi-reahstic two- 
nucleon S3 potential of Afnan and Tang jl^, as modified in Ref. In the present work we neglect 

the three-nucleon potential. For the A-nucleon (AN) potential we have chosen that proposed by Bodmer 
and Usmani jT9|. We have investigated the effects of many-body forces involving the A by explicitly 
considering a ANN interaction or, alternatively, by including a density dependence in the bare AN force. 

The interactions generate strong NN and AN correlations, beyond those described by a mean field 
approach. In our work the correlation effects are accounted for by a Jastrow correlated wave function, 
where the two-body correlations depend on the interparticle distances only. More realistic hamiltonians 
demand for a strong spin and isospin dependence in the correlation factor. However, we consider the 
simple Jastrow ansatz as a first step towards a more complete description of the hypernuclei. The same 
strategy has been followed for doubly closed shell nuclei, where the accuracy of the variational method 
has now reached the same level as in nuclear matter [|o| . 

The FHNC equations for a A embedded in infinite nuclear matter were developed in Ref. . In Ref. 
p2t , doubly closed shell nuclei in Zs-coupling were studied by FHNC and Jastrow-type correlations, while 
the J j-coupling was introduced in such a scheme in Ref. P3| . It was so possible to study nuclei ranging 
from ^^C to ■^"^Pb within the correlated basis function theory, using cluster expansion at all orders and 
avoiding such shortcuts as low order cluster truncations and local density approximations. The results 
obtained in those papers are the basis of the present FHNC study of A-hypernuclei. 

The plan of the paper is as follows: section |l| is devoted to a brief description of the hamiltonian and 



of the correlated wave function; section |in| gives an outline of the FHNC theory for the hypernucleus; in 
section 



. [V the expressions for the variational energy are presented and discussed; the results are given in 
section ^and the conclusions are drawn in section Vl . 



II. HAMILTONIAN AND TRIAL WAVEFUNCTIONS. 

The system under study is a hypernucleus composed by A nucleons and one A hyperon. The formalism 
we have developed is however general enough to describe any nuclear systems containing a non-nucleonic 
particle, treated as an impurity in the nucleonic fluid. 

We write the hamiltonian of the system as: 



Hn+a = Hn + -ffA = 



2mi 



i—l j>i—l 

where we have considered two- and three-body interactions between the hyperon and the nucleons and 
only two-body ones between the nucleons. In the hamiltonian (0) the purely nucleonic part, Hf^, has 
been separated from that involving the A, H^^. The separation will be useful in the remainder of the 
paper. 

As already stated in the introduction, we use the modified NN S3 potential, which is fixed to reproduce 
the s wave scattering nucleon-nucleon phase shifts at low energies, gives reasonable results both for light 
nuclei and nuclear matter and allows for utilizing only central two-body correlations between the nucleons. 
The operatorial structure of S3 is: 

4 

V''''{^,J) = Y.VPin,)0':^ , (2) 

where 0^~^~^ = 1, <Ti-<Tj, Ti- Tj, <Ti ■ (TjTi - Tj. Modern nucleon-nucleon potentials are more complicated 
and contain the essential tensor terms, absent in S3. On the other hand, we are mainly interested in 
presenting and testing the formalism. Adopting the S3 potential is sufficient for this purpose and has the 
advantage of avoiding the technical complications deriving from the use of more realistic potentials. 

Within the same philosophy, we have chosen the A-nuclcon potential proposed by Bodmer and Usmani 
p9[. This potential has central, spin and space-exchange components: 
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2) = Vl^{r,2) + Vl^{r,2)fTK ■ cr^ + V!^^{r,2)P. , (3) 

with 

Vlj,{r)=v^{r){l~e) , (4) 

VUr)^\{v,-v,)n{r) , (5) 

VX^M=e«oM, (6) 
where Px is the spatial exchange operator and uq is defined as: 



"oM = . , -K, - 7("^ + i^t)n{r) . (7) 



l + exp(i^) 4^ 



The one-pion exchange term is 



TAr) =(l + - + j^) ^ (l - , (8) 

\ fir [firy J fir \ / 

with the pion mass /i=0.7 fm^^ and the cutoff parameter c=2.0 fm~^. The values of the various parame- 
ters VsjVt, e, Wc, R, a are given in Ref. [^4| . In this reference, a three-body A-nucleon-nucleon interaction 
was added to the just described A-nuclcon one in order to fit the hypernucleus ^J^O empirical Ba by the 
variational Monte Carlo (VMC) method. The ANN interaction has an important, attractive two-pion 
exchange part, which is the source of a tensor term. The tensor contribution, however, results to be zero 
if simple central correlations are adopted. As a consequence, such a semi-microscopic ANN interaction 
is too repulsive in a Jastrow correlated model, as we shall discuss later. For this reason, we simulate 
the effects of three- and more-body potentials by using an approach similar to that developed by Fried- 
man and Pandharipande for pure nucleonic matter |2^]. We multiply V^^{1, 2) by a density dependent 
function 

2) = exp{-7[p(ri) + p(r2)]/2} , (9) 

where p{r) is the nucleonic density. The value of the parameter 7 is fixed to reproduce the empirical 
value of a single A binding energy in nuclear matter, B^'^ = 30 MeV. 

The binding energy of the A particle is defined as minus the difference between the energies of the 
nuclear systems with and without the A : 



Ba — „ . , -, < Hn+a >a+a — < Hn >a 
2.7 + 1 ^ 



^ ^ ' < r^+A >A+A -<n^>A] , (10) 



2i + 1 ^ 



where we have used the notation: 



< X >y= J drnX^y , (11) 

with Y = A + A 01 A. In the previous equations we have averaged on the third component of the A total 
angular momentum and we have explicitly singled out the difference between the center of mass energies 
of the two systems. 

Our approach is based upon the variational principle, therefore the energy functionals given in eq. (p^ , 
< Hjfj^A >A+A and < >a, should be independently minimized with respect to variations of the 
many-body wavefunctions, ^'y. The hypernucleus wavefunction is chosen as: 



\^A+A >= (jlfAirM)^ ^^i^„,irA)\^A > , (12) 

where 0^;^-^ is the mean-field single particle wavefunction of the A in jj-coupling, fAi^Ai) is a two-body 
scalar (Jastrow) correlation between the hyperon and a single nucleon and > is the correlated wave 
function describing the remaining A nucleons. In our calculation this function is defined as: 
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\^A>^ 1^ n fN(^v)j <fN^z{l,...,A) , (13) 

where ^n,z is the Slater determinant ol a set of single particle wave functions of N neutrons and Z 
protons. Obviously, A = Z + N. 



III. DENSITIES AND FHNC EQUATIONS 

The energy expectation values are evaluated by means of the Fermi hypernetted chain cluster 
diagrams resummation technique. As already mentioned, our approach consists in considering the A as 
an impurity in the nucleonic fluid. The FHNC equations for an impurity in homogeneous matter were 
derived in Ref. pl[| for a A hyperon in symmetric nuclear matter and in Refs. [EGp^] for a '*He atomic 
impurity in liquid "^He. 

Our aim is the evaluation the A binding energy in finite nuclear systems. To this purpose, we extend 
the FHNC equations developed in Ref. for doubly closed shell nuclei in the j j-coupling scheme with 
a Jastrow correlated wave functions. 

The basic quantities to be analyzed are the one- and two-body densities (OBD and TED). It is conve- 
nient to define densities for each type of particle, protons p, neutrons n and A: 

A 

PlA+Air) = (E -^(r - rfe)n">A+A (14) 

k=l 

p^(r) = (5(r-rA))A+A (15) 

where a = p,n, P" is the projector operator over the particle of a type. An averaged sum on the third 
components of the A particle angular momentum is understood. As in Refs. |p^ ^ p3|| the densities are 
evaluated by using cluster expansion techniques. The nucleonic OBD are divided in two parts: 

PlA+Air) = pUir) + Pt Air)- (16) 

The first part is given by the sum of the cluster diagrams containing only nucleons, and represents the 
nuclear bulk contribution to the OBD. One of the diagrams contributing to this part is the A diagram of 
Fig.^. The second part of eq. called rearrangement term, is obtained by summing all the diagrams 
containing the A as an internal particle, like the diagram B of Fig.^ and provides the modification of the 
nucleon OBD due to the presence of the A impurity. A mean-field description of the hypernucleus does 
not provide any rearrangement term. Finally, the diagrams where the A is an external particle, like the 
C diagram of Figj^, are summed in the A-OBD, p^(r). The separation of the one-body densities in bulk, 
rearrangement and A terms will be very useful for the calculation of the A binding energy in the nucleus 
and it will be exploited also for the calculation of other quantities occurring in the FHNC scheme. 
In analogy with the one-body densities, we define the two-body densities as: 

A 

pt.A+Air, r') - ( E -^(^ - - rOO|,P,"Pf )^+A (17) 

= P2,lA(i",r')+P2iA(r,r') , 

where a, (3 = p,n and is one of the operators defined in eq.(^. The AN two-body densities are 
defined as: 

A 

p^,^(r,r') = {S{v-r^)Y,S{r'-rk)Ol,P^)A+A , (18) 
fe=i 

A 

P^%iv,v') = {5{v-v^)Y^5{v' -vu)P^,P^)a+a , (19) 
fe=i 

where = 1, cta • cr^, arc the operators of the AN interaction (||), and Px is the space exchange 

operator between the A and the k-nucleon. 
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In the FHNC approach the TBD are described in terms of two-body correlation functions fN{fki) and 
fAi^Ak), and of uncorrelated densities. The latter are built from a single particle basis, defined in the 
present paper as a set of jj coupled wave functions of the type: 

= ^nljmir) = R'nljir) J2 < ^^^^^'"^ > : (20) 

where t = p,n,A. Yi^^ are the spherical harmonics, < l^^s\jni > are the Clebsch-Gordan coefficients 
and Xs the two-component Pauli spinors. 

In the nucleonic case, the uncorrelated two-body densities are: 

p[?(ri,r2)=^</.r(ri)</.^(r2) 

fe 

= ^pr'"(ri,r2)xl(l)x.'(2). (21) 

s.s' 

The explicit expressions of the parallel, s — s' , and antiparallel s = —s' densities are given in Ref. p3[ 
together with the other purely nucleonic FHNC quantities. 

We briefly discuss here the quantities involving the A. The uncorrelated A one-body density is given 
by: 

Po(r) = ^E'^"y™(r)0il,^.„(r) = ^ (i?,^,,.(r))' . (22) 

The one-body A uncorrelated density matrix is needed to calculate diagrams where the spatial coordi- 
nates between the A and one nucleon are exchanged. Its expression is: 

p^(ri,r2) = ^R^i^{n)R^,i^{r2)Pi{coseu) , (23) 

where we have indicated with Pi the Legendre polynomial of Ith degree and the angle between the 
two vectors. 

A detailed description of the FHNC equations in the nucleonic matter, in a jj coupled single particle 
basis computational scheme, has been done in Ref. We do not give here the equations for this case, 
but we rather present and discuss the changes related to the presence of the A. The FHNC A-densities 
are: 

Piir)^C^ir)^eir)p^ir) , (24) 
e(0 = ^cxp[L/,^(r)] , (25) 

Ca - 1 dr p^ir) expiU^ir)] , (26) 

p^j(ri,r2) = C^in) (e,^(r2)5^i'(ri,r2) +Cf (r2)ffi''(ri,r2)) , (27) 
p^|(ri,r2)-0. (28) 

The functions C^e(?'), U^{r) and g^^{ri,r2) are analogous to their nucleonic partners whose expressions 
are given in the Appendices A and B of Ref. |2^, and they are obtained with the substitution a ^ A. In 
extending the expressions given in the reference, we should remember that there is no sum on A, since it 
is just a single external impurity, and that we consider the same correlation function between A and any 
type of nucleon: fAf3{ri2) = /a(7'12) with l3 ^p,n. 

Eq.(p8[) shows that, in our case, the spin component of the two-body density, P2 2' zero. The reason 
lies in the fact that only diagrams corresponding to exchanges between identical particles contribute to 
the spin component of the two-body density for Jastrow correlated wave functions. As a consequence 
they do not affect the AN-TBD. Direct diagrams, where the particles are not exchanged, may contribute 
to P22 only if the correlation contains spin dependent components (which is not our case). Moreover, 
this type of contribution has been found to be generally small p4| . 
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The evaluation of the diagrams containing the space exchange AN potential requires the introduction 
of a space exchange density, given by: 

P2,x(r,r') = ^/l(|r-r'|)p^(r,r')p^^"(r,r') , (29) 

where p**"(r,r') is defined as in Appendix B of |23), with the substitutions fai3{'ri2) ~^ /A(''i2)/iv(»'i2) 
in the construction of the gi^d functions. 

The expressions of the rearrangement terms are quite lengthy since the A may be present in any of the 
elements involved in the densities. We give below the expression of the contributions to the one-body 
density: 

P?.aW = (UZAir) [P^ir) + U:{r)] + U^^^)) exp[;7,"(r)] (30) 
= airWSAr) + CirWZAr) = ^aM , (31) 
CAr) = U2Ar) eMU2{r)] ■ (32) 

The U^^j^ g functions are given in Ref. |^^, while the C/^ and the rearrangement parts of the two-body 
densities are given in the Appendix. 

IV. THE VARIATIONAL ENERGY 

In this section we discuss the evaluation of the variational energy of a A hypernucleus in the framework 
of the FHNC approach. The potential energy is separated into its nucleonic and A pieces: 

< V;v+A >A+A=< Vn >a+a + < ^a >a+a • (33) 
The expectation values can be expressed in terms of the two-body distribution functions as: 

4 



<Vn>a+a = \ ^ J2 J c'i"idr2F«(ri2)p2iA+A(i"i,i-2) 
= <Vn>a+1 E E/'^'-i*2T^'(^i2)P2iA(ri,r2) , (34) 
<Va>a+a = E E/ rfridr2FAV(ri2)p^,^(ri,r2) + 

a.—p,7i q—1 

^ f dr,dr,V^Ari2)pti^i,r2) , (35) 

a=p,n 

where < V > indicates the average on the third components of the A angular momentum and < Vn >a 
is the bulk nucleonic potential energy. 

Because the minimizations of < Hn+a >a+a and < Hn >a must be, in principle, carried on inde- 
pendently, the nucleon-nucleon Jastrow correlation, fNir), might be different in the two cases, as well as 
the nucleon single particle wave functions, 0fc(r). However, we found that in our cases /Af(?') and 0fc(r) 
do not practically change in going from the A to the A -\- A system. As a consequence, in our approach 
the nucleonic part of the hypernucleus wave function has been kept the same as in nucleus. This fact 
allows for obtaining the potential energy contribution to the binding energy of the A by subtracting the 
pure nucleus < Vn >a from < Vn >a+a- The remaining part is separated in two contributions: the 
potential energy due to the interaction A-nucleon, or interaction energy, and the modifications of the 
nucleon-nucleon potential energy due to the presence of the A, or rearrangement energy. We put in 
evidence these contributions by rewriting: 

Va =< Vn+a >a+a -<Vn >a= Vl + Vf , (36) 

with 
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a— p,n Vg— 1 



Aa I 



and 



The separation in interaction and rearrangement terms can also be done for the A kinetic 
evaluated here by means of the Jackson-Feenberg expression p2| . Following an analysis similar 
done for the potential energy, we obtain: 

7a =< Tn+a >a+a - <Tn >a~ tI + . 
The interaction kinetic energy is given by: 



^A 



4rnA 



drdp^i(r)e(r) 



with 



1 

47r 



i?™;,!'^!) A%('^l)- 



^a+l)oA 



The expressions of t[f\] and -D^y can be obtained from Ref. [g3| (Eq.(16) and Appendix C). 
The rearrangement part of the kinetic energy is: 



rriR 

^A 



^ E / C^I"lC^I"2t[/w(ri2)]p2a,A(l"l'J'2) 

a.f3—p.n 

^T. f rfrip?i(ri)e,A(ri) 

|1 ^ y'dridr2p?2(ri,r2)[C(ri)C(r2)5?d"A(ri,r2)+ 



(Ce,A(ri)C(r2) + C"(ri)C"A(r2)) gZir^r^)] 



2m 



/ dridra p^^g p(ri, r2)i/"^_p_A(ri, r2) 



where 



+PT3,A(ri,r2)i7"c,A,A(''l''^2) 



c%,A(ri,r2)) =4"a('^i)[C(^2) ((5r(ri,r2) - 1) 7V,",,^(ri, r2)) + ^^(^^^^(ri, rs) 
+ {C{r2) - l)A^i^5(ri,r2) +4"(r-2)5r(ri,r2R",^Z5(ri,r2y 
+ C(ri)ee"A('^2)5r(ri,r2) (7Vrc,D(ri,r2) + i?,",,c(ri, rs)) 
+ C(^i)[Ce"(^2) ((5^(1-1, r2) - l)iV,",^,,,A(ri,r2) 
+ ^^^(ri, r2)iV,",,c(ri, r2) + iVi^^^^Cn, rs)) 
+ (C(r2)-l)ivi^,'^,A(ri,r2) 

+ C ('^2 ) A (ri , r2 )i?e",D (i-i , r2 )+ gr (ri , r2 ) i??,,z,, A (i-i , r2 ) ) 



(37) 



(38) 



energy, 
to that 



(39) 



(40) 



(41) 



(42) 



(43) 
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Again the expressions of p^^, for k = 1,2,3 can be found in the appendix C of 
The difference of the center of mass kinetic energies is given by: 

rpcm ™A ^ mcm ^ " 

J A — ^ 



Am + fUA 

and the expression of < T™' >a is given in eq. (23) of Ref. te3| 



< ^^"^ >^ 'JTT^ V / *iA(ri) , (44) 

4(Ato + rriA) J 



V. RESULTS 



We have studied the ground state structure of those A-hypernuclei having N and Z values corre- 
sponding to doubly closed shell nuclei in jj-couphng, ^^C, ^JO, ^^^a, ^^Ca, and ^"^Pb. The 
same set of nuclei where the neutron in the highest neutron single particle level has been substituted 
by a A hyperon {^J^C, ^^O, '^^Ca, ^^Ca, ^°Zr and ^"^Pb) have been also considered. In order to deal 
in FHNC with the partially occupied neutronic level, we adopt the following procedure: the degeneracy 
factor 2j + 1, multiplying the partially occupied state contribution to the uncorrelated nucleonic OBD, 
has been replaced by the factor 2j. In this way, the OBD results normalized to A — 1. Moreover, in 
the spin parallel part of the uncorrelated nucleonic TBD, , the 2j + 1 factor has been substituted by 
2a/ {21 + 1){21 + 1/2) — 4Z + 2j — 1, while the (much smaller) antiparallel part, p^, remains unchanged. 
This choice ensures a correct normalization of the TBD. 

Two different mean field potentials have been used to build the single-particle nuclear wave functions. 
The first potential is a harmonic oscillator (HO) well with the same constant, b^, for protons and 
neutrons; the second choice is the Woods-Saxon (WS) potential used in the calculations of Ref. |Q. The 
parameters of the WS potential are different for protons and neutrons. We have always used a harmonic 
oscillator potential for the A single-particle potential. 

The correlation functions have a gaussian form: 

fx=AM{r) ^l~ax exp(-/3xr2) . (45) 

In order to calculate Ba, we first minimize the bulk energy of the nucleus, < H]\[ >a- The minimization 
has been performed for each nucleus in two ways, depending on the nucleonic single-particle used. In 
the case of the HO potential, the nuclear correlation, /Ar(r), has been taken from nuclear matter and the 
nucleus energy has been obtained by varying only the values of the oscillator constant of the HO mean 
field. For the Woods-Saxon case, the mean field potential has been kept fixed as in Ref. and the 
parameters of the correlation have been varied. We consider the WS model as the most realistic one, since 
its parameters have been determined to reproduce at best the nuclear densities. It must be noticed that, 
in the region of the variational space corresponding to the energy minimum, the energy itself is rather 
insensitive to small changes of the WS parameters. Therefore it is possible to reproduce the densities 
and the radii without excessively spoiling the quality of the minimum. 

For the description of the nucleonic part of the hypernucleus we take the same /7v('") and mean field 
potential as in the corresponding nucleus. Then, < H^i+a >a+a is minimized by changing only the 
A-HO constant and the parameters of the A-nucleon correlation function, fA{r)- 

The parameters of the nuclear correlation ( ^ ) with the HO mean field are aN — 0.7 and [3n — 2.0 fm^^. 
The nuclear matter binding energy per nucleon with this correlation and the S3 potential, at saturation 
density {pnm = 0.16 fm^'^), is i3jvM=14.43 MeV, close enough to the empirical value, -BArA/,cmp=16 MeV, 
and comparable with the best, more sophisticated potentials on the market. For each nucleus the value of 
the oscillator constant, ^at, has been fixed to get the energy minimum. In Table ^ we compare our HO and 
WS binding energies and rms charge radii with their experimental values. A general reasonable agreement 
is found, considering the relative simplicity of both wave functions and interactions. In particular, the 
WS radii are close the experimental ones. The ^^C nucleus represents somehow an anomaly, having a 
good estimate of the radius but a very small binding energy. The origin of this disagreement is still under 
investigation. It may be related to the inadequacy of a spherical model description of this nucleus. 

The parameter 7 of the density dependent potential (DDP), described in Section II, has been chosen 
to reproduce the A binding energy in nuclear matter, Ba,nm- We find 7 — 2.2 br? and Ba,nm — 30.26 
MeV and this value of the parameter 7 has been used in the remaining finite systems calculations. The 
parameters of the A-HO mean fields and of the /a(?') correlations at the < Hn^a >a+a minimum are 
shown in Table ||. 
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In Table Til we compare the contributions to the A binding energy in the doubly closed shell hypernuclei 
calculated with the HO nucleonic mean field with those evaluated in nuclear matter We have also separated 
the interaction (I) and rearrangement (R) terms in the kinetic and potential energies. The upper part 
of the table gives the results with the two-body AN interaction only and the related A binding energy, 

. In the lower part of the table the contribution of the three-body ANN interaction of Rcf. pi] and 

(3) 

the corresponding binding energy J5^ are shown. Finally we present the results obtained with the DDP 
potential. 

Since the empirical value of B\ in is estimated to be about 13 MeV the result obtained with 

the AN force only is too attractive. The VMC calculation of Ref. [H gives S^^|^mc("O)=27-5±2.0 
MeV, in qualitative agreement with our value. The empirical binding energy was reproduced in the VMC 
approach by the inclusion of an explicit ANN force. Using the same three-body potential, we obtain 
b'"^\^\0)=&.81 MeV. We have already pointed out that central, Jastrow, correlations underestimate the 

attractive contribution of the ANN interaction, as clearly appears from the B\ value. This potential, 
when used in conjunction with Jastrow correlated wave functions, does not even bind the A in heavy 
nuclei. Tensor-like correlations are needed in order to effectively use potentials induced by one- and 
more-pion exchanges in spherically symmetric systems . The introduction of a DDP, fitted to the A 
binding in nuclear matter, brings i?A reasonably close to the empirical estimate, even with a relatively 
simple correlation. All the results presented hereafter have been obtained with the DDP. 

In Table IV we give the binding energies of a A in its Is ground state, calculated with the HO single 
particle potential for several hypernuclei and nuclear matter. We explicitly show the overall interaction 
and rearrangement contributions. The difference of the center of mass kinetic energies has been included 
in the rearrangement part. The increase of the A binding energy along A is mostly produced by the 
interaction energy and, specifically, by its potential energy part. In contrast, the dependence of the 
rearrangement energy on A is much weaker. We found an analogous behavior the A energies in the Ip 
and Id states. To complete the information we give in Table the results for the A energies in the Ip 
and Id states with the nucleonic HO mean field. In Table [Vl| we compare the A binding energies obtained 
with a Woods-Saxon nucleonic mean field with the experimental energies. For the ^^^r the comparison 
is done with the energies measured in *^Y. 

The results for the A binding energies are summarized in Fig .B, where they are presented as a function 
of A'"^/^ and compared with the experimental energies of Refs. ||] (dots), @1 (triangles) and [H for the 
^^C. The figure gives the A energies in the Is, Ip and Id states. The quality of the agreement with the 
experiment is rather good, in spite of the simplicity of the model. By inspecting the figure in more details, 
we find that all the calculations underbind the A in Carbon by ~ 30%, and provide a steeper variation of 
_Ba with respect to the experiments. This behavior resembles that already described in Carbon nucleus. 
It is also worth noticing that, as expected, the heavier nuclei, Zr and Pb, are better described by the 
Woods-Saxon well than by the HO one. 

The effects of the AN correlations on the one-body densities are shown in Fig.|, where the proton 
densities in four nuclei are compared. The dotted lines are the IPM densities, the dashed lines are the 
densities obtained in a purely nucleonic FHNC calculations, and, finally, the full curves represent the 
proton densities obtained when a A hyperon is added in the s wave to the doubly magic nucleonic core. 
In general, correlation effects are more important in Oxygen and Calcium than in the other two nuclei. 
Moreover, the presence of the A does not heavily modify the nucleonic densities. This is better shown in 
Fig.^, where the differences between proton (panel a) and neutron (panel h) densities with and without 
the A are given. The differences shown in the figure have been amplified by a factor 1000. The heavier 
nuclei seem to be much more stable against deformations produced by the presence of the A. 

This picture might change when the A is inserted in the nucleonic core by substituting the last bound 
neutron. The differences, multiplied by 1000, between the nucleonic densities of the isotopic hypernuclei 
with A-l-1 and A hadrons are given in Fig|5| The panels show the differences between the protons (upper) 
and neutron (lower) densities. The conservation of the charge and mass numbers for each hypernucleus 
implies that the curves in the upper panel are normalized to zero, while those in the lower panel are 
normalized to one. 

The analysis of this figure indicates that the nucleonic part is more perturbed if the A substitutes a 
neutron having a low angular momentum, like in O, Ca, and Pb. In Zr the neutron modified into a A is 
lying on the 139/2 level (^=5), and its wave function is peaked at the surface of the nucleus, as shown in 
panel (&) of Fig.H. In any case, heavy nuclei have again densities more rigid against deformations induced 
by the hyperon. 
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The A densities in ^\0, ^\Ca, ^aZt and ^'^^Ph for s and p waves are shown in FigJ|. The dashed hnes 
are the IPM densities and the full ones those obtained in the full calculations. Results very similar have 
been obtained for the ^%0, "^^Ca., ^^Zi and ^^APb hypernuclei. Again, the differences between the IPM 
and FHNC approaches are smaller for the heavier systems. 

VI. CONCLUSIONS 

In this work we have studied some properties of hypernuclei containing one A hyperon in the framework 
of the correlated basis functions theory. The A binding energies in the s, p and d states, its density and 
the rearrangement part of the one-nucleon densities have been computed for hypernuclei whose nucleonic 
core has a doubly closed shell structure in jj coupling. The hyperon has been either added to the core 
(i^C, iJO, ^iCa, 4^Ca, ^\Zt and ^ospb) or substituted to a neutron in its highest shell (i^C, ^^O, "^Ca, 
48Ca, 90Zr and ^ospb). 

The correlated wave function contains gaussian, Jastrow correlations for the NN and AN pairs. The 
correlations act on a hypernuclear shell model wave function generated i) by a harmonic oscillator or 
Woods-Saxon mean field for the nucleonic part and ii) by a harmonic oscillator well for the A single 
particle potential. Cluster expansion and Fermi hypernetted chain resummation technique have been 
used. The energy of the system has been minimized with respect to variations on the parameters of the 
wave function. We have employed semirealistic hamiltonians with the S3 two-nucleon potential of Afnan 
and Tang, the Bodmer and Usmani AN potential and either a model of ANN three-particle interaction, 
still proposed by Bodmer and Usmani, or a density dependent modification of the AN force, in order to 
take into account many-hadron interactions. 

Using the ANN potential with Jastrow correlations severely underestimates the A binding energy, since 
the important, attractive tensor component of the two-pion exchange potential does not contribute in 
this model. In contrast, the density dependent modification of the AN potential, fitted to the A binding 
in nuclear matter, provides results for the A binding energy in nuclei in encouraging agreement with the 
experimental data. In fact, the disagreement is less than 10% in all hypernuclei, for both s and p A states. 
The only exception is provided by the carbon hypernuclei, showing a discrepancy with the experiments of 
^^30% that could be ascribed to the inefficiency of a spherically symmetric description of these systems. 
Moreover, our analysis stresses the importance of adopting a good independent particle wave function as 
a starting point. In fact, the Woods-Saxon nucleonic mean field, fitted to the experimental one-nucleon 
densities, gives a better description of the A binding than the harmonic oscillator model, especially in 
the heavy hypernuclei. 

The nucleonic core polarization effects due to the presence of the hyperon are more important in the 
lighter hypernuclei than in the heavier ones, both for the energy and for the one-nucleon densities. The 
rearrangement contributions to the A binding energy go from ~ 20% in C to ~ 5% in Pb and nuclear 
matter. As far as the A-density is concerned, the influence of the correlations is much more visible in O 
and Ca than in the heavier Zr and Pb. 

The results shown in this paper have been obtained by means of relatively simple hamiltonians, which 
entitle to use equally simple wave functions. However, hyperon-nucleon interactions may now be built on 
more microscopic grounds. It is to be expected that an effective use of these hamiltonians will ask for a 
more sophisticated correlation, with a strong dependence on the relative state of the correlated pair. The 
situation closely resembles that in doubly closed shell nuclei and in nuclear matter, where the introduction 
of modern potentials has prompted the extension of the correlated basis functions and FHNC theories in 
that direction. Another interesting field of application of this technology is the study of hypernuclei with 
two or more hyperons, especially in view of a better determination of the hyperon-hyperon interaction. 
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APPENDIX 



In this Appendix we collect the different contributions to the rearrangement part of the one- and 
two body densities, where one of the internal particles must be a A hyperon. Following this rule, we 
obtain: 

USA^i) = I rfr2|^^(r2) W(ri,r2) - ^?j^(ri,r2) - 5,«j^(ri, r2)(5?i^(ri, r^) - 1)) 

+ Ef^^(^2)(x,tA(ri,r2) - <A(ri,r2) - 5,tA(ri, r2)(5f (n, r^) - 1) 

-5,t(ri,r2)5,tA(ri,r2)) (46) 
+ef (^2)(x«f^(ri,r2) - <A(ri,r2) - SXj,{v,,v,){g^^{r,,r,) - 1) 
-^f (n, r2)<A(ri, r2) - S^,^^^{rur,)gl^ {ri,r,) 
-5,t(ri,r2)5,tA(ri,r2)) 

+CiA(r2)(x,t(ri,r2) - i?,t(ri,r2) - S^,^ {r,,r,)ig"J {r,,r2) 1)) 
<A(r2)(^,f (ri,r2) - <(ri,r2) - S^,^ {v„v,){g^^ {v„v,) - 1) 

-^.t(ri,r2)5f(ri,r2))j| , 

t^e':A(n) = / rfr2|^^(ri)(x«/(ri,r2) - i;e"/(ri,r2) - 5«j^(ri, r2)(5?j^(ri, r2) - 1) 

-S2^{r„r,)g:^\r„r,)) (47) 
+ Ef^^^(^2)(^etA('-i''-2) - <A(ri,r2) - 5,t^(ri,r2)(5f (ri,r2) - 1) 

-5^^(1-1, r2)ffrf|A(ri,r2) - S^l^{ri,r2)g"^{ri,r2) 
-5d,f(ri,r2)5e|A(i"i,r2)) 
+ef (r2)(x:i^A(ri,r2) - <A(ri,r2) - <A(ri, r2)(5,t (ri, r2) - 1) 
-5,f (ri,r2)5,tA(ri,r2) - 5etA(ri, ^2)5^ (ri, r2) 
-^et(ri,r2)5,tA(ri,r2) - <A(ri,r2)5:i'(ri,r2) 
(ri,r2)5etA(ri,r2) - 5,tA(ri, r2)5ef (ri, r2) 
-5,t(ri,r2)5etA(ri,r2)) 
+ef,A(^2)(x,f (ri,r2) - i?:^ri,r2) - S:,^ {r„r2){g"J {r^r^) - 1) 
-^."^(ri,r2)5et(ri,r2)) 

+ef,A(^2) (n, r2) - (ri, r2) - (n, r2)(5,t (ri, r2) - 1) 
-^et(ri, r2)5f (ri, r2) - (n, r2)5:,^(ri, r2) 



-5,t(ri,r2)5ef (ri,r2) 



+ 4C(»'2)(5'^c,p(ri,r2)5"c,p,A(ri,r2) + _p_A(ri, r2)5f^c_p(ri, r2) 



+ 
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'S'"c,A(ri, r2)5cc,A,A(ri, r2) + -S^cA.ACri, r2)5cc,A(ri, r2)) 
- 2(p?_p(ri,r2)Ar,V,A(ri,r2) +p^,^(ri,r2)7V«,^,A(ri,r2)) 

+ 4C,A(r-2)(5e",p(ri,r2)5c",p(ri,r2) + 5cV(ri> r2)5cV(ri> ^2)) | • 

The rearrangement parts of the two-body distribution functions are: 

P2,lA(ri,r2) = ed"(ri)(5dtA(ri,r2)C,^(r2) (ri,r2)C,V(r2) + 

5rftA(ri, r2)^f (r2) + fff (ri, r,)^^,^) 
+ e?,A(r-i)(fff (ri,r2)ef (r2) +5f (ri,r2Xf (r2)) 
+ C(ri)(ffetA(ri,r2)?f (r2) +5:i'(ri,r2)eiA(r2) + 

5etA(ri,r2Kf(r2)+.g,f (ri,r2)^f,A(r2)) 
+ C,A(r-i) (5:<f (ri, r2)$f (r2) + fl^f (ri, r2)ef (r2)) , 
<„,A(ri,r2) = A^l,A(ri,r2) +^l,A(ri,r2) , 
5l,A(ri,r2) = A^l,A(ri,r2) +X«^„_^(ri,r2) , 
S.tAC'-i' r2) = 5.".^(ri, r2)<^(ri, r2) , 

5.tA(ri, r2) = 5,",^(ri, r2)<^(ri, r2) + (ri, r2)<^(ri, r2) , 
5etA(ri,r2) = 5,",^(ri,r2)<A(ri,r2) +5f (ri,r2)B,tA(ri,r2) 
+ 5:<f (ri, r2)B,tA(ri, r2) + (n, r2)S,tA(ri, r2) 
- 25a/3(5cc,p(ri,r2)Bcc,p,A(ri,r2) +5cc,A(ri>r2)-B^c,A,A(ri,r2 
5cc,i5,A(ri, r2) = 5r(ri, r2)B« ,^,A(ri, r2) + g^,^n{ri,r2)B^^,A{ri,r2) • 
The equations for the nodal diagrams are: 

^l,A(ri,r2) = (X;^^(ri,r3)ed^r3)|xf (r3,r2) + iVf (r3,r2)) + 
(^:iA(ri,r3)a'(r-3)|^,t(r3,r2) +7V,Tf (r3,r2)) + 
(x«](ri,r3)e,],,A(r3)|X,Tf (r3,r2) +iV,Tf (r3,r2)) + 

with m,n = d,e and 

^'(^3) = ^dd.H'eKrs) + (1 - SM,w)a{r3) , 
Qwi'^s) = SddM'Q^Ai^s) + (1 - Sdd,ii')Qj^{r3) . 

Finally, the equations for the cc nodal are: 

^cc,D,A(ri,r2) = (^c",D,A(ri,r3)ee(r-3)|.9S=,p(r3,r2)) + 
(X,",,^(ri, r3)C,A(^3)|5c" ,p(r3, r2)) + 
(X«_^(ri, r3)C(r-3)|5c"c,p,A(r3, r2)) + 



EE 
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-N. 



{P)c 



■D,A 



:,£;(ri,i'3)CA('^3)|.gcc,A(i"3,r2; 



■SD (^c",£;,A(i"i,r3)C(^3)l5cc,A(r3,r2) 

(x,",,s(ri,r3)C(^3)|<?cc.A,A(r3,r2) 
(ri,r2) = (pS(ri,r3)C"A(''3)l5?c,p(r3,r2)) + 

(pS (ri , r3)C (^^3) l^e" ,P,A(r3, r2) + N^SA^^^^^)'^ 
(pS(ri,r3)(C(r3) - l)l<^;,A(r3,r2)) + 
SD (pE(ri,r3)C"A(r3)|5°^^(r3,r2)) + 

(p^(ri,r3)C(r3)|^e%,A(r3,r2) + iv(:!4"A(r3,r2)) + 
(p^(ri,r3)(C(r3)~l)|<^XA(r3,r2)) ^ , 
where P,A. li D = P, then S = A and sp = -1; if £) = A, then E = P and = 1. 
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HO 


ws 


exp 




biv 


R 


B/A 


R 


B/A 


R 


B/A 




1.56 


2.48 


2.63 


2.47 


2.78 


2.47 


7.680 


16 Q 


1.48 


2.43 


7.07 


2.69 


6.29 


2.73 


7.976 


«Ca 


1.66 


3.08 


9.00 


3.29 


8.12 


3.48 


8.551 




1.62 


3.00 


7.57 


3.35 


6.79 


3.48 


8.666 


»Zr 


1.74 


3.57 


10.07 


4.09 


7.30 


4.04 


8.710 


208 pj^ 


2.05 


4.67 


10.24 


5.52 


8.03 


5.50 


7.867 



TABLE I. Optimum values of the nucleonic oscillator constant, fojv (in fm), of the rms radii R (in fm) and of 
the binding energy per particle B/A (in MeV) for different nuclei. The Woods-Saxon potential parameters are 
those of Ref. [||. 





HO 


WS 




bA 




/3a 


bA 


OA 


/3a 




1.92 


0.80 


3.0 


1.92 


0.80 


3.1 




1.86 


0.80 


2.9 


1.86 


0.80 


2.9 


^Ca 


2.00 


0.75 


2.9 


2.04 


0.75 


2.9 


^Ca 


2.00 


0.75 


2.9 


2.04 


0.75 


3.0 


^IZr 


2.24 


0.70 


2.9 


2.48 


0.80 


3.1 


20lPb 


2.68 


0.70 


3.0 


2.96 


0.80 


3.1 



TABLE II. Optimum values of the A oscillator constant, 6a (in fm), and of the parameters of fA(r), qa and 
/3a (in fm^'^) for the two nuclear mean fields used in this work. In nuclear matter, aA=0.70 and /3a =3. 2 fm~^. 



Hypernucleus 




''iCa 


^A^a 




^"«Pb 


NM 


Ti 


17.12 


18.98 


20.78 


19.96 


19.24 


10.13 


rpR 

A 


-1.66 


-2.65 


-3.63 


-3.87 


-3.52 


-0.49 


rj-'Cm 

A 


0.46 


0.14 


0.13 


0.06 


0.02 




vl 


-42.76 


-61.63 


-69.78 


-80.54 


-87.87 


-59.89 




2.31 


3.59 


4.82 


5.08 


5.09 


2.22 




24.53 


41.57 


47.68 


59.31 


67.04 


48.03 




17.72 


37.80 


49.12 


74.41 


77.72 




4) 


6.81 


3.77 


-1.44 


-15.10 


-10.68 






12.75 


21.81 


26.99 


34.65 


37.89 


17.77 


Br" 


11.78 


19.76 


20.69 


24.66 


29.15 


30.26 



TABLE III. Contributions, in MeV, to the A binding energy for various hypernuclei calculated within the HO 
model and compared with the nuclear matter (NM) results. See text. 
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Si 


Ba 


Ba 


a'^ 


9.54 


-1 


85 


7.69 


13 


10.16 


-1 


74 


8.42 




12.41 


-1 


39 


11.02 




12.89 


-1 


11 


11.78 


^Ca 


20.76 


-1 


18 


19.58 


^Ca 


20.84 


-1 


08 


19.76 


^Ca 


21.90 


-1 


39 


20.51 


^iCa 


22.01 


-1 


32 


20.69 


^XZr 


25.88 


-1 


31 


24.57 


'XZr 


25.93 


-1 


27 


24.66 




30.72 


-1 


60 


29.12 




30.74 


-1 


59 


29.15 


NM 


31.99 


-1 


73 


30.26 



TABLE IV. Interaction, rearrangement and total binding energies, in MeV, of a A in the Is state with the 
HO nucleonic mean field. Also the energy of a single A in nuclear matter is shown. 





Ba{Ip) 


BA(ld) 




0.73 






1.38 




^Ca 


10.01 


0.85 


*lCa 


10.32 


1.13 


JCa 


11.08 


2.00 


^Ca 


11.27 


2.19 


«XZr 


16.85 


8.69 


«iZr 


16.95 


8.80 


20|Pb 


23.54 


17.18 


20lPb 


23.58 


17.24 



TABLE V. A binding energies, in MeV, in the Ip and Id state calculated with HO nucleonic mean field. 







Ba(1s) 




Ba(1p) 




BA(ld) 




th 


exp 


th 


exp 


th 


exp 


A^ 


7.57 


10.8±0.1 




0.1±0.5 






13 n 


8.30 


11.7±0.1 




0.8±0.5 






lo 


11.21 


12.5±0.4 


1.12 


2.5±0.4 








11.99 




1.76 








^XCa 


19.61 


18.7±1.1 


10.27 


10.1±0.3 


1.37 


1.0±0.5 


^Ca 


19.96 




10.61 




1.64 




^Ca 


21.14 




11.89 




2.98 




ICa 


21.34 




12.10 




3.19 




^XZr 


23.19 


22.1±1.6* 


16.76 


16.0±1.0* 


9.96 


9.5±1.0* 


«lZr 


23.27 




16.87 




10.08 




20|Pb 


27.52 


26.5±0.5 


22.78 


21.3±0.7 


17.33 


16.5±0.5 


20lPb 


27.55 




22.82 




17.39 





TABLE VI. -Ba, in MeV, for a A in the Is and Ip states calculated with the WS nucleonic mean fields. The 
experimental binding energies are from Refs. [p9|-|3l|. The energies labelled with an asterisk have been measured 
for the nucleus. 
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a 



A 



a 




a 




a 




A 



B 



C 



FIG. 1. Examples of cluster diagrams contributing to the one body densities. The A diagram contributes 
to p".A(r), the B diagram to pi^/^ir) and the C diagram to Piir). The dashed lines represent the dynamical 
correlations, fo{rij) — 1, {D = N in K and D = A the rest of the cases) and the oriented lines the statistical 
correlations, po{ri,rj). A black dot associated with a point implies integration over its coordinates. The label a 
indicates the nucleon. 
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FIG. 2. A binding energies for the Is, Ip and Id states, as function of A~^^^. The white circles are the energies 
calculated using harmonic oscillator (upp er p anel) and Woods-Saxon mean field potentials (lower panel). The 
experimental energies are taken from Ref. |29[ (dots) and (triangles). The full lines connecting the theoretical 
values have been drawn to guide the eyes. 
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FIG. 4. Differences between proton (panel a) and neutron (panel h) densities with and without core polariza- 
tion. 
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FIG. 5. Differences between proton (panel a) and neutron (panel b) densities of isotopic hypernuclei with A+1 
and A hadrons. 
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